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IJ  ABSTRACT 


Thir  contract  has  been  concerned  vith  the  application  of  ultrashort  optical 
puls  to  high  resolution  imaging  optical  radar.  Previous  work  has  demonstrated 
two  ! in ens ional  range-doppler  imaging  of  a  laboratory  target.  This  was  accom¬ 
plished  by  range  gating  the  returned  signal  with  a  short  duration  local  oscillator 
pulse.  This  technique  is  capable  of  extremely  high  time  resolution,  but  it  suf¬ 
fers  from  the  problem  of  inefficient  use  of  the  received  signal  energy.  An  alter¬ 
nate  signal  processing  technique  involving  frequency  domain  sampling  has  been 
ilucussed  briefly  in  previous  semiannual  report.  During  this  period  further 
generalization  of  this  scheme  has  been  carried  out.  This  scheme  is  considerably 
more  versatile  than  the  initial  description  indicated;  it  is  capable  of  performing 
several  very  useful  signal  processing  operations.  It  can  perform  a  time  scaling 
on  repetitive  ultrafast  optical  signal  and  stretch  its  time  duration  to  any 
desired  duration  while  retaining  the  functional  form  of  the  signal.  In  addition, 
it  can  synthesize  an  optical  match  filter  for  an  arbitrary  signal,  loth  of  these 
operations  can  be  performed  with  electronics  that  are  much  slower  than  the  time 
scale  of  the  original  signal. 


This  report  reviews  breifly  the  application  of  ultrashort  pulses  to  imaging 
radars,  and  shows,  tne  need  for  an  alternate  signal  processing  scheme.  A  math- 
ematical  description  and  physical  interpretation  of  frequency  domain  sampling  is 
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presented  and  its  application  to  time  scaling  and 
matched  filtering  is  discussed.  An  experiment  to 
demonstrate  the  technique  is  outlined  and  will  be 
carried  out  during  the  next  reporting  period. 
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TECHNICAL  REPORT  SUMMARY 


This  contract  has  been  concerned  with  the  application  of  ultrashort  optical 
pulses  to  high  resolution  imaging  optical  radar.  Previous  work  has  demonstrated 
two  dimensional  range-doppler  imaging  of  a  labo.'ntory  target.  This  was  accomplished 
by  range  gating  the  returned  signal  with  a  short  duration  local  oscillator  pulse. 

This  technique  is  capable  of  extremely  high  time  resolution,  but  it  suffers  from  the 
problem  ol  inefficient  use  of  the  received  signal  energy.  An  alternate  signal  pro¬ 
cessing  technique  involving  frequency  domain  sampling  has  been  discussed  briefly  in 
the  previous  semiannual  report.  During  this  period  further  generalization  of  this 
scheme  has  been  carried  out.  This  scheme  is  considerably  more  versatile  than  the 
initial  description  indicated;  it  is  capable  of  performing  several  very  useful 
signal  processing  operations.  It  can  perform  a  time  scaling  on  repetitive  ultrafast 
optical  signal  and  stretch  its  time  duration  to  any  desired  duration  while  retaining 
the  functional  form  of  the  signal.  In  addition,  it  can  synthesize  an  optical  matched 
filter  for  an  arbitrary  signal.  Both  of  these  operations  can  be  performed  with  elec¬ 
tronics  that  are  much  slower  than  the  time  scale  of  the  original  signal. 

This  report  reviews  briefly  the  application  of  ultrashort  pulses  to  imaging 
radars,  and  shows  the  need  for  an  alternate  signal  processing  scheme.  A  mathematical 
description  and  physical  interpretation  of  frequency  domain  sampling  is  presented 
and  its  application  to  time  scaling  and  matched  filtering  is  discussed.  An  experi¬ 
ment  to  demonstrate  +  he  technique  is  outlined  and  will  be  carried  out  during  the 
next  reporting  period. 
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SECTION  1 

IMAGING  OPTICAL  RADARS 


This  contract  has  been  concerned  with  the  application  of  extremely  short  dura¬ 
tion  optical  pulses  to  imaging  radars.  The  problem  that  has  been  considered  is  that 
of  imaging  a  target  whose  size  and  distance  is  such  that,  it  cannot  be  resolved  by 
conventional  optical  imaging.  If  there  is  relative  motion  between  the  target  and 
the  radar  transmitter,  range-doppler  processing  of  the  received  signal  may  be  used 
to  generate  a  two  dimensional  image  of  the  target.  The  situation  is  illustrated 
schematically  in  Figure  la.  T  le  target  is  taken  to  be  a  collection  of  scattering 
centers  that  are  rotating  as  a  rigid  body  about  an  axis  of  rotation.  If  the  received 
energy  is  analyzed  in  t^me  and  frequency,  a  range-doppler  image  of  the  target  may  be 
constructed  as  shown  in  Figure  lb.  The  contours  of  constant  range  are  lines  perpen¬ 
dicular  to  the  line-of-s ight  between  the  transmitter  and  the  target,  while  the  con¬ 
tours  of  constant  doppler  shift  are  parallel  to  the  line  of  sight.  Together,  they 
impose  a  coordinate  system  on  the  target  which  is  used  to  generate  -he  image.  In 
this  simplified  case,  the  image  looks  exactly  like  the  target.  For  more  complicated, 
three  dimensional  targets,  the  situation  is  more  complicated,  but  the  range-doppler 
image  still  gives  the  projection  of  the  target  on  the  down-range  and  cross -range 
directions . 

The  solution  in  the  down-range  direction  is  determined  by  the  time  resolution 
of  the  radar;  i.e.,  by  the  reciprocal  oi  the  bandwidth  of  the  transmitted  waveform. 

If  the  bandwidth  is  B,  then 


Ax 


c 

2B 


The  resolution  in  the  cross  range  direction  depends  on  how  well  the  doppler  frequencies 
can  bn  l-" solved.  The  doppler  shift  for  a  scatterer  located  at  position  y  is 

_  2^ Qv 
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where  U  is  the  rotation  rate  and  u>  is  the  optical  frequency.  If  the  target  is  observed 
for  a  time  t,  then  the  minimum  resolvable  dop;  ler  shift  Is  2rr/t.  We  have  then 

2  a;  1  Ay  >  2n 
t 

2  TT  c  _  _X _ 

>  2u)Dt  ’  2A0 


R M 


f 


■■  »■  .m  i  m. 


M-920UT9-1+2 


where  9  =  Qt  is  the  angle  through  which  the  target  has  turned  in  the  time  T.  Both 
Ax  and  Ay  are  independent  of  the  range  of  the  target  in  contrast  to  the  resoa  tier 
of  a  conventional  imaging  system  where  the  minimum  resolvable  element  size  an  rea^e. 
linearly  ^ ith  range.  The  magnitude  of  the  doppler  shift  is 

U>d  =  Urryft 

\ 

For  a  target  size  of  1  meter 

U3 ,  =  Itt  x  10^  ft  at  X  =  l.o6  p, 
d 

=  In  x  10 at  \  =  10.6(1. 

Laser  bandwidths  of  ~  1  GHZ  are  presently  available  in  COg  lasers  at  10. 6u,  and 
it  is  expected  that  much  larger  bandwidths  will  be  available  in  the  future.  Neodymium- 
YAG  lasers  have  bandwidths  that  could  give  time  resolutions  of  ~  >0  psec.  A  time 
resolution  of  100  psec  corresponds  to  a  range  resolution  of  1.5  cm,  and  1  Psec  to 
0.015  cm.  Optical  radars  are  thus  potentially  capable  of  extremely  high  spatial 
resolution,  In  order  to  take  advantage  of  this  very  high  resolution,  means  for  pro¬ 
cessing  extremely  fast  optical  transients  must  be  developed.  Most  of  the  'fort  on 
this  contract  has  been  concerned  with  this  problem.  Although  the  bandwidth  of  the 
transmitted  signal  is  necessarily  very  large  if  high  resolution  is  desired,  the  actual 
information  rate  is  not,  The  information  desired  from  the  received  signal  is  an  image 
of  the  target.  If  we  assume  that  we  want  an  image  with  100  resolvable  elements,  each 
with  a  10  bit  gray  scale  and  that  we  obtain  the  image  in  1  millisecond,  then  the  actual 
information  rate  is  ~  106  blts/sec.  This  is  far  less  than  the  bandwidth  of  the  trans¬ 
mitted  waveform  which  could  be  10-^  Hz  or  higher .  This  observation  leads  to  the  en¬ 
d-union  that  with  proper  signal  processing,  the  desired  image  could  be  obtained  with 
a  relatively  slow  detection  system. 

A  variety  of  waveforms  could  be  used  to  obtain  an  image  of  the  target.  As 
mentioned  above,  the  iown  range  resolution  is  determined  by  the  reciprocal  ol  the 
bandwidth  B  of  the  signal  and  the  cross  range  resolution  by  the  total  duration  r 
of  the  signal.  The  "average  resolution"  c\n  be  taken  as 

(AxA  y)2  ~  (BTT 

so  that  all  signals  with  the  same  time-bandwidth  product  are  equivalent  a"  far  as 
resolution.  One  possibility  is  to  use  a  periodic  train  of  short  pulses  that  can  be 
obtained  conveniently  from  a  mode-locked  laser.  This  situation  is  illustrated  in 
Figure  2a.  The  time  between  the  pulses,  if  derived  from  a  mode-locked  laser,  would 
typically  be  a  few  nanoseconds.  The  received  signal  will  consist  ol  a  quas i-periodi 
train  oi  longer  duration  pulses  as  shown  in  Figure  2b.  Each  of  these  pulses  is  a 
narrow  band  signal  in  times  of  the  optical  center  frequency,  and  can  be  described  in 
terms  of  an  envelope  f (t)  and  an  instantaneous  phase  *(t).  Adjacent  received  pulses 
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are  nearly  identical  but  there  is  a  slow  variation  of  f  and  $  that  occurs  over  a 
time  corresponding  to  the  reciprocal  of  the  doppler  spread.  The  corresponding  fre¬ 
quency  domain  description  of  the  transmitted  and  received  signals  is  shown  in  Fig¬ 
ure  3. 


One  way  of  processing  this  signal  is  to  range  gate  it  and  measure  the  amplitude 
and  phase  within  each  range  element.  Frequency  analysis  of  the  amplitude  and  phase 
within  each  range  element  gives  the  doppler  spectrum  within  that  range  element.  Mea¬ 
surements  on  each  range  element  may  then  be  used  to  constr'.  ct  the  range-r’oppler 
image  of  the  target.  The  signal  could  be  processed  by  the  scheme  shown  ii  Figure  4a. 
The  signal  is  heterodyned  against  a  local  oscillator,  detected,  and  range  gated.  If 
the  local  oscillator  is  at  the  same  center  frequency  as  the  signal,  a  quadrature 
channel  with  the  local  oscillator  shifted  by  \/4  is  also  needed.  If  there  is  a  fre¬ 
quency  offset,  this  is  not  necessary.  For  a  stationary  target,  the  output  of  each 
range  element  will  vary  at  the  offset  frequency  and  will  be  proportional  to  the  sig¬ 
nal  received  in  that  range  element.  If  the  target  is  rotating,  the  doppler  spectrum 
from  each  range  element  will  be  superimposed  on  the  signal.  In  this  system,  the 
detectors  and  the  range  gating  electronics  must  be  sufficiently  fast  to  accommodate 
the  desired  time  resolution  of  the  received  signal. 

An  alternate  scheme  is  shown  in  Figure  4b.  Here  the  signal  is  heterodyned 
against  a  mode-locked  local  oscillator  which  may  be  derived  from  a  frequency  shifted 
version  of  the  transmitted  signal.  Heterodyning  only  occurs  when  the  pulsed  local 
oscillator  is  present  so  that  the  local  oscillator  provides  the  necessary  range 
gating.  Different  optical  propagation  paths  may  be  used  to  shift  the  position  of 
the  range  gate  as  shown.  This  system  does  not  require  fast  detectors.  The  time 
resolution  is  determined  by  the  duration  of  the  local  oscillator  pulse.  The  detectors 
need  only  respond  to  the  if  frequency. 

An  experiment  was  carried  out  to  demonstrate  this  type  of  imaging.  The  experi¬ 
ment  is  shown  schematically  in  Figure  5.  The  source  used  was  a  iw  mode  locked 
neodymium-YAG  laser  that  emitted  pulses  of  ~  100  psec.  duration  and  with  a  separation 
of  —  3  nsec.  The  beam  was  expanded  to  a  diameter  of  ~  6  cm.,  and  illuminated  the 
target.  The  polarized  component  of  the  return  signal  was  reflected  from  the  Glan- 
Tnompson  prism  and  was  shifted  in  frequency  by  4o  MHz  by  an  acoustic  modulator.  It 
was  then  heterodyned  against  a  local  oscillator  derived  from  the  transmitted  beam 
as  shown. 

The  target  consisted  of  a  bar  holding  two  retroref lectors  that  were  separated 
by  a  distance  of  3.5  cm.  The  target  was  rotated  at  0.5  rps  about  a  horizontal  axis. 
The  sw^ep  of  the  spectrum  analyzer  was  initiated  when  th°  bar  was  at  45°  with  the 
upper  half  receding  from  the  transmitter,  and  was  completed  before  the  target  had 
rotated  appreciably.  The  target  was  slowly  scanned  in  range,  and  a  spectrum  was 
taken  at  corresponding  points  in  each  revolution.  This  slow  scan  in  range  allowed 
the  various  range  elements.  of  the  signal  to  be  measured  sequentially  and  eliminated 
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the  need  for  the  parallel  channels.  The  resulting  signal  was  displayed  on  a  storage 
scope.  The  horizontal  position  of  the  beam  vas  scanned  in  synchronism  with  the 
spectrum  analyzer  sweep  and  the  vertical  position  was  scanned  along  with  the  target 
motion.  The  output  of  the  spectrum  analyzer  modulated  the  beam  Intensity.  This 
arrangement  produced  a  range-doppler  Image  directly.  Such  an  image  Is  shown  in 
Figure  6.  The  upper  portion  shows  the  results  with  the  laser  operating  on  a  single 
mode  (no  range  resolution)  and  the  lower  portion  shows  the  results  with  the  laser 
mode-locked.  The  excellent  spatial  resolution  (~  2  cm)  Is  evident. 


Although  this  type  of  signal  processing  is  capable  of  excellent  time  resolu  1  -  , 
It  is  inefficient  in  the  use  of  the  received  signal  energy.  Even  for  the  case  shown 
in  Figure  4b,  where  parallel  range  samples  are  taken,  much  of  the  signal  is  was  e 
due  to  the  fact  that  it  must  be  spilt  among  the  separate  channels  for  each  range 
element.  This  Is  a  fundamental  problem  In  any  time  sampling  system.  Efficient  use 
of  the  signal  would  require  some  kind  of  a  fast  switch  or  beam  deflector  to  switch^ 
the  various  range  elements  of  the  signal  to  different  detec  ors  or  '° 
a  fast  detector  capable  of  resolving  the  time  variation  of  the  received  s  gnal  As 
has  been  mentioned  above,  the  actual  Information  rate  of  the  signal  is  relatively 
low  so  there  should  exist  techniques  to  process  the  signals  that  do  not  require  id 
bandwidths.  One  method  Is  sampling  In  the  frequency  domain  rather  than  In  the  1 
domain  ^hls  type  of  processing  is  discussed  In  detail  in  the  next  section.  I 
“  „lL  £»  the  fact  that  It  Is  possible  to  spatially  separate  the  frequency 
components  of  an  optical  signal  vlth  a  passive  device,  a  spectron,vter.  This  vlll 
allow  a  measurement  of  the  Fourier  Transform  of  the  optical  signal  In  an  e 
way .  It  Is  not  possible  to  perform  the  same  separation  of  the  time  segtnen  s  o 
signal  without  a  fast,  active  element.  The  frequency  domain  sampling  also  offer* 
the  possibility  of  performing  additional  signal  processing  operations. 
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FREQUENCY  DOMAIN  SAMPLING 


2.1  Arnuyn  In 


We  will  first  consider  the  application  of  frequency  domain  sampling  to 
the  measurement  of  a  repetitive  train  of  extremely  last  optical  signals.  The 
signal  can  be  represented  as  a  Fourier  series 


n  =  N 

S(t)  =  23  sn  cos  [(UJ0  +  nO  )t  + 

n  =  N 

s(t)  cos  (uj0t  +$(t))  -T/2Z0ZT/2 

2  rr 

repeats  every  T  =  -rr- 


2.1 


The  quantities  that  are  to  be  measured  are  s(t),  the  instantaneous 
amplitude,  and  $(t),  the  phase  of  the  signal.  It  will  be  assumed  throughout 
that  the  signal  is  narrow  band  in  terms  of  the  optical  center  frequency  u>0, 
although  it  may  be  very  wide  band  in  terms  of  detector  capabilities.  The 
Fourier  coefficients  are  given  by 


T/2 

Sn  sin  <pn  =  -l/T  f  sft)  sin  (nQt  -  $(t))  dt 
-T/2  ' 

T/2 

Sn  cos  cpn  =  +l/T  |  s(t)  cos  /nQt  -  $(t)j  dt 
-T/2 


2.2 


The  general  form  of  the  signal  processor  is  shown  in  Figure  7.  The 
signal  is  first  combined  on  a  beam  splitter  with  a  reference  waveform.  This 
waveform  will  be  denoted  by 


n  =  N 


R(t)  =  r^  cos  |(u)q  +uj^  +  nQ)t+i>n 


n  =  N 


yr(t )  cos  (uj0  +  u>L)  t  +  0(t)|  -T/2ZQZT/2 
I  repeats  every  T  =  2  n / Q 


2.3 
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This  reference  waveform  is  periodic  with  the  same  period  as  the  original 
signal,  but  its  optical  certer  frequency  is  shifted  by  an  amount  u>]_.  It  is 
assumed  that  tu^ZZ  nZZu)0*  The  composite  signal  leaving  the  beam  splitter  i, 


n  =  +N 


C(t)  =  (Sn  C°S  l^°  +  n°‘  ^  +  *«]  + 

rn  cos  |(“o  +  «>i  +  n  0  )  t  +  )]) 


2.h 


n  =  -N 


The  spectrum  of  this  signal  is  shown  schematically  in  Figure  7-  This  signal 
is  then  passed  through  a  bank  of  narrow  band  optical  filters.  These  filters 
are  centered  on  the  lines  of  the  spectrum  of  C(t).  (The  lines  are  actually 
doublets  as  shown  in  Figure  7  ,  but  the  spacing  of  the  doublet  is  small  by 
assumption).  The  filters  will  all  be  assumed  to  be  identical.  The  transmission 
of  the  m  th  filter  for  the  n  th  spectral  line  will  be  denoted  by  Tm  (%  )  e  'vrr‘>n/, 
Since  the  filters  are  identical,  the  transmission  of  a  given  line  depends  only 
on  its  displacement  from  the  center  frequency  of  the  filter  so  that  the  trans¬ 
mission  may  be  written 


T  (m-n)  e 

The  output  of  the  m  th  filter  is  given  by 


i  •if  (m-n) 


Fm  =  £  T(m-n)  |Sn  cos  [(t»0  +  nQ  )  t  +  cpn  +  t  (n-m)] 
n 

+  r  n  cos  [(u>0  +  U)]_  +  nfl)  t  +i?n  +  i|r  (n-m)] 

Ea^h  of  these  outputs  is  the  square  law  detected  to  give  a  signal 

Gm  r  £  £  T(n-m)  T  (n'-m)  jsn  cos  [(tl)0  +  nfl)  t  +tpn  +  i|i(n-m)| 
n  n' 


2.5 


Is 


n 

t 

n 

t 

n 


cos 

cos 

cos 


|(<l)0  +u>]_  +  nfi)t+#n  +  ♦  (n-m)]  |  . 
[(«0  +  n'  0)  t  +  spn,  +  t(n'  -m)] 

[( U>0  +  u»i  +  n'  n)  t  -  i>n.  +  ♦  (n'  -m)]| 


2.6 


The  detectors  do  not  respond  at  the  optical  frequency  so  that  only  the  difference 
terms  in  the  products  of  cosines  are  significant.  This  gives 
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Gm  =  n  S  T(n-m)  T(n'-m)  /  SnSn*  cos  J  (n-n'  )  fit  +  cpn-cpn '  +  i|l(n-m)  -*(n-m)| 
n  n'  ' 

+  rnrni  cos  [(n-n')flt  +  t?n  -  t>ni  +  ty(n-m)  -  ^(n'-m)] 

2.7 

+  Snrni  cos  |(n-n' )  Ot  -  +  cpn  -  t?n «  +  lit(n'-m)  -  ty(n-m)] 

+  Snirn  cos  |(n-n')Ot  +  uj]_t  -  cpn  1  +  i?ni  +  ty(n-m)  -  ^  ( n  *  -  m  )| 

Each  of  these  signals  is  now  filtered  with  a  narrowband  electrical  filter 
centered  at  u)^.  This  eliminates  all  terms  in  the  above  expression  except 
those  for  which  11 '  t  n.  The  resulting  signal  is 

Hm  “  E  T(n-m)?  Snrn  c°s  K*  ’  *n  +  *n)  2.8 

n 

Here  the  phase  factor  of  the  optical  filters  has  dropped  out.  Each  of  these 
outputs  is  now  multiplied  by  cos  (tug  +  mflg  where  u>g  an^  ^2  are  chosen  at 
will  (but  with  f^ZZuig).  They  are  th-n  summed  to  give 


w(t) 


5  EE  T(n-m)P  Snrn  {  cos  |(ti^  +  u)2  +  Mfig)t  -  cpn  +  ^n] 
n  m 


+  cos 


u)^  +  M  fig  ^ 


2.9 


This  may  be  written  in  a  more  convenient  form  by  letting  p  =  n-m  and  eliminating 
m.  The  result  is 

W(t)  =  \  (E  T(P)"  cok  P°2 t)(E  Snrn  cos  [(®1  +  *2  +  n^2  )t  -  «>n  +  !’n]) 

((<«1  +  u>2  +  )t -cpn  +  *n]) 

K  -'«1  +  nfi2)t  +  Vn-^n]) 
[(u)2  -  +  nflgjt  +  cpn  -  djj 


\  (E^p)2  sin  pn2t)(E  Snrn  sir 
'  p  ' '  n 

\  (Et2(P)C  cos  pfig^E  Snrn  cc 

\  (ET^)?  sln  Pn2t)(E  Snrnsi 


We  now  consider  some  special  cases.  If  the  filters  are  sufficiently  narrow 
that  they  only  pass  one  spectral  line,  then  T(p)  =  6(p)  and 


-8- 


M-92DUT9-1+2 


W 


00  =  \  ESnrn  cos  [(<“1  +  <*>2  +  nn2)  t  '  +  *nl 

+  \  EVn  cos  |(“l  +  «2  +  n02  )  t  +  cpn  -  i?J 


2.11 


Suppose  now  that  the  reference  waveform  consists  of  a  train  of  extremely  short 
pulses.  In  this  case  rn  is  essentially  constant  and  &n  =  0.  We  have  ther 

W(t)~  |  £Sn  cos  [(uj1  +  uj2  +  nn2)  t  -  cpj 
+  |  ZSn  cos  [(u^  -  uj2  +  n02)  t  +  cpj 


2.12 


Reference  to  Eq.  2«1  will  show  that 

1 


W 


(t)  =  -  S(-at)  cos  [(oi1  +  u>2)  t  .  i  (-a  t )] 
+  |  S(of  t)  cos  |(u)1  -  u)2)  t  +  $  (at)] 


2.13 


=  w+  (t)  +  w_  (t) 
for  -T2/2Z0ZT1/2 


repeating  with  period  T,,  ■  211  /»12 
Q2 

and  with  a  =  — 

The  second  term  is  this  an  exact  replica  of  the  original  optical  signal  but 
with  a  time  scale  that  has  been  stretched  by  a  factor  Q/iV>.  The  center  fre¬ 
quency  of  the  signal  has  also  been  changed  from  the  optical  frequency  m0  to 
an  r.f.  frequency  u)^  -  (l^.  The  first  term  is  a  time  reversed  version  of  the 
optical  signal,  stretched  in  time,  and  with  a  center  frequency  U)-j_  +  U)2.  The 
signal  processor  has  extracted  the  amplitude  and  phase  of  each  of  the  compo- 
nets  of  the  line  spectrum  of  the  original  optical  signal  and  has  imposed  them 
on  a  line  spectrum  at  an  r.f.  frequency.  This  synthesizes  the  original  signal. 
The  signal  and  its  time  reversed  version  may  be  readily  separated  by  the  proper 
choice  of  u)^  and  U)2  together  with  frequency  filtering. 

bet  us  now  consider  the  more  general  case  where  rn  and  j?n  are  arbitrary. 
The  component  at  ujg  -  is 

W(t)  =|  £Snrn  cos  [(m1  -  u>2  +  nQ2  )  t  +  tpn  -  2.14 
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The  cross  correlation  of  the  signals  s(t)  and  r(t)  is 
T/2 

J  ^  S(t)R(t-T)dt  B  2  n  Snrn  COR  l^o  +  n°)  T  +  <Pn  -  fln]  2.15 

The  f inction  W_(t)  may  be  se^n  to  be  a  repetitive  time  display  of  the  cross 
•or relatloc  ’unction.  The  function  W+(t)  is  just  a  time  reversed  version  of 
this  '■lgnai.  In  the  case  where  S(t)  and  R(t)  are  identical,  the  processor 
>-TIv.«,,  the  autocorrelation  function  as  a  function  of  time;  i.e.,  it  behaves 
a  1  matched  filter.  In  the  case  where  R(t )  is  a  train  of  pulses  of  finite 
iuration  pulses  rather  than  extremely  narrow  pulses,  the  original  signal  is 
tv<  0 1 ve  1  with  a  degraded  time  re  olutlon  given  by  the  duration  of  the  pulses 
■  >  1  tltutiug  R(t). 

We  may  now  return  to  the  more  general  case  given  in  Eq.  2.10.  Let 
T(p  Y'  -os  pfl2t  =• 

Eo 

T(p)  sin  pQ2t 

We  note  that  ii  the  bandpass  of  the  optical  filters  are  symmetric;  i.e.,  if 
|T(P)|  -  |Tfp)|,  then  Ts(t)  =  0. 

G(t)  =  E  snrn  cos  [(uj1  +  u>2  +  nfl2  )  t  -  cpn  +  t?n] 

=  g(t)  cos  j(un  +  ajj)  t  +  y  (t)]  2.17 

(periodic  at  T0  =  LVr/^) 

lor  n  narrowband  waveform,  it  will  be  true  that 


T0(t) 

Ts(t> 


Z,  snrn  sin  [(*1)^  +  U)2  +  nO  )  t  -  q>n  +  i?n]  =  g(t)  sin  |(u»x  +  u)2)  t 

+  V(t )]  2'V 

So  that  we  have  for  the  uj  +  ui-  co-iponent  of  Eq.  2.10 


W+  =  T,(t)  g(t)  cos  [(u^  +  tt,)  t  +y(t)J  +  Tc(t)g(t)  sin  [(m  +  w? )  t 

+  V  (t)l 

I  2.19 
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A  similar  expression  holds  for  W_.  The  shape  of  the  filter  transmission  func¬ 
tion  multiplies  the  output  of  the  system  by  a  known  function  of  time.  The 
result  that  would  be  obtained  with  very  narrow  filters  can  be  computed  from  the 
result  obtained  from  wider  filters.  Alternately,  the  filter  functions  could  be 
chosen  to  perform  additional  signal  processing  operations. 

The  entire  discussion  above  has  been  concerned  with  repetitive  signals. 

The  frequency  domain  sampling  technique  is  not  restricted  to  this  case.  It 
may  be  applied  to  the  measurement  of  a  single  optical  transient.  The  processing 
is  somewhat  more  complicated  in  this  case  because  it  corresponds  to  heterodyning 
with  zero  difference  frequency  and  care  must  be  taken  to  avoid  spectral  over¬ 
lapping.  A  system  to  perform  this  measurement  is  illustrated  schematically  in 
Figure  9-  We  assume  again  a  signal  and  a  reference 

R(t)  ‘  jh)i  I  '1“‘t  r<<«>  * 

The  signals  are  combined  on  a  beamsplitter  to  produce  a  composite  spectrum 
s(u>)  =  s(u>)  +  r(u>) 

The  two  signs  are  obtained  from  opposite  sides  of  the  beamsplitter.  The  com¬ 
posite  signal  is  then  passed  through  a  filter  bank.  At  the  output  of  each 
filter  a  signal 

F(»)  =  R  (S(w)  +  r(*))  Tm(u>) 

is  obtained.  This  is  square  law  detected  and  integrated  over  time  to  obtain 


Gm  =  ls(®)  -  r^)i2  lTrn(U})|2  du> 

=  r.{l8M2±*e  (r(uu)s* (uj))  +  |r(u))|2j  Tn(»)  dcu 


The  difference  between  the  plus  and  minus  signals  is  then  taken  to  give 


Um  =  2  Sh  Re  (r(®)s*(®))  Tm(U}^  du} 


-VI- 


M-9?oU79-^2 


If  the  filters  are  extremely  narrow,  Tn(tu)  as  6  (a)-a>m )  and 

U  =  2  Re(r(u)  )  s*(ti)  )) 

m  v  m  m  < 

If  the  same  operation  is  repeated  with  the  reference  wave  shifted  by  X/4,  the 
output 

u„  ■  2  T" (>■(*,>  *  *(»,)) 

•will  be  obtained.  We  thus  obtain  a  sampled  version  of  r(uj)  C  (u>).  This 
quantity  is  the  Fourier  Transform  of  the  cross  correlation  function;  .e. 

f  R(t)  S  (t-T)  dt  =  n  eiU)t  r(iu)  s  (u>)  du> 

u  —08  o  _0 

These  parallel  outputs  may  be  used  to  reconstruct  the  signal  at  any  desired 
frequency  anl  repetition  rate  in  a  manner  identical  to  the  previous  lscussion. 
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2.2  Discus5 ion  am  Experimental  Implementation 


Tie  analysis  of  the  previous  section  described  the  operation  of  the  fre¬ 
quency  sampler  in  the  time  domain.  Figure  8  shows  the  situation  in  the  fre¬ 
quency  domain,  and  illustrates  why  the  turn  frequency  domain  sampling  has  been 
used.  The  various  components  of  the  line  spectra  of  R(t)  pass  through  the 
filters,  are  isolated,  and  are  distributed  to  separate  detectors.  If  an 
appropriate  local  oscillator  were  available  at  each  detector,  the  amplitude 
and  phase  of  each  of  the  components  could  be  measured.  This  would  completely 
characterize  the  signal.  Just  such  a  set  of  local  oscillators  can  be  obtained 
by  filtering  a  periodic  train  of  short  pulses  through  the  same  filter  bank  as 
the  signal.  Each  channel  can  be  operated  in  the  quantum  nrise  limited  regime, 
and  all  of  the  energy  receive]  in  each  channel  is  used. 

An  alternate  way  of  looking  at  the  operation  ->f  the  system  is  to  consider 
the  response  of  one  of  the  narrow  band  filters  to  a  short  optical  pulse.  The 
output  of  the  filter  has  a  duration  corresponding  to  the  reciprocal  of  the 
bandwidth  of  the  filter.  Thi3  should  be  made  comparable  to  the  duration  of 
the  signal  that  is  to  be  measured.  The  reference  signal  can  thus  heterodyne 
with  the  signal  over  its  entire  duration.  The  precise  details  «>"  the  result  ng 
signal  are  given  above. 

We  may  now  consider  means  of  implementing  such  a  detection  system.  One 
of  the  first  considerations  is  the  resolving  power  required  of  the  filters. 

Let  us  assume  that  the  line  shape  of  the  filter  is 


+  II .  The  modulating  functions  are 


This  is  a  porlodic  function  with  period  T2  =  2 tt  /Cl2 *  We  may  rewrite  this  in 
the  form 

p  "  P  2 

T  (t)  =  1  +  2  y  sln^ttp  cos  p02t 
O’2  ,  P2. 

p  =  1 


'(p> 


sin  op 
op 


where  the  first  zeros  occur  at  op  = 


T  (t)  =  0 


p  =  I’ 


T0(t)  -  £ 

p  --P 
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For  values  of  a  <tt  this  simply  represents  a  triangular  pulse  extending  from 
-  (y/0  <  t  <a/0  as  shown  in  Figure  10.  When  y  =  tt  the  pulse  extends  over 
thp  entire  interval  -  T/2  <  t  <  T2/2.  For  larger  values  of  a,  the  base  line 
lifts  as  shown  in  the  figure.  For  the  particular  values  a  =  RN  ,  the  function 
Tc(t)  =  1.  The  reason  for  this  is  shown  in  Figure  10.  Although  the  filter 
has  xinite  transmission  off  its  center  frequency,  its  transmission  for  the 
lines  of  the  spectrum  of  S(t)  is  zero  except  for  the  line  on  center. 

For  ease  in  data  interpretation,  the  filters  should  be  made  quite  narrow 
(or  they  should  satisfy  the  condition  illustrated  in  Figure  11).  They  should 
not,  however,  be  made  so  narrow  that  they  are  capable  of  resolving  the  if  fre¬ 
quency,  since  this  will  pervent  the  heterodyning  of  the  signal  ani  the  refer¬ 
ence  waveforms.  The  bandwidth  Auj  of  the  filters  should  be  in  the  range 
uJif  <Acu<  0  .  For  a  typical  mode  locked  laser,  0/2n~  3  x  Hz,  at  a 
wavelength  of  1.06  p,,  U30/2tt  =  3  x  10^  Hz.  This  corresponds  to  a  resolving 
power  of  u>0/Au)  ~  10^.  This  resolving  power  is  easily  obtainable  with  a  Fabev- 
Perot  etalon.  It  can  be  obtained  with  more  difficulty  with  a  grating,  a  large 
grating,  possibly  with  several  passes  would  be  required.  The  grating  approach, 
however,  has  several  advantages  which  will  be  discussed  later. 

The  number  of  filters  required  is  determined  by  the  desired  time  resolu¬ 
tion  (which  is  ultimately  determined  by  the  bandwidth  of  the  reference  signal). 
Figure  12  shows  the  response  of  the  system  to  an  amplitude  step  for  a  filter 
bank  employing  5>  11  and  21  filters.  Figure  13  shows  the  amplitude  and  phase 
of  the  response  of  the  system  to  a  pulse  of  constant  amplitude  and  linearly 
increasing  phase  (a  frequency  offset)  for  a  filter  bank  with  11  components, 
and  Figure  lU  for  one  with  21  components.  Figure  15  shows  the  response  to  a 
chirped  pulse  for  11  and  21  components.  These  responses  were  computed  by  cal¬ 
culating  the  Fourier  coefficients  of  the  applied  waveform,  truncating  them 
at  the  desired  number  of  components  and  then  reconstructing  the  series  with 
the  truncated  set  of  coefficients. 

The  problems  involved  in  the  experimental  demonstration  of  the  frequency 
domain  sampling  techni  iue  are: 

1.  attainment  of  sufficient  resolving  power  in  the  optical  filter  bank 

2.  simultaneous  detection  of  the  signals  from  each  of  the  filters 

3.  modulation  of  the  output  of  each  detector  by  the  appropriate 
modulation  waveform 

recombination  of  the  modulated  outputs  to  regenerate  the  time- 
scaled  version  of  the  original  signal. 
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The  first  of  these  requirements  determines  the  type  of  spectrometer 
that  is  required.  The  spectrometer  must  have: 

1.  a  resolving  power  of  the  order  of  10 ^ 

2.  high  throughput 

3.  produce  an  output  of  all  the  frequency  components  of  the  incident 
beam  simultaneously. 

A  Fabry  Perot  etalon  can  satisfy  the  first  two  of  these  requirements  readily. 
It  does  not,  however,  satisfy  the  third.  When  a  Fabry  Perot  etalon  is  operated 
with  a  collimated  input  beam,  it  acts  as  an  optical  band  pass;  i.e.,  it  passes 
the  wavelength  to  which  it  tuned,  with  high  transmission,  and  it  reflects  all 
the  other  wavelengths  with  low  loss.  If  it  is  operated  with  an  input  beam  from 
an  extended  source  or  one  with  more  divergence  than  the  diffraction  limit  of 
its  aperture,  then  it  acts  as  a  spectrometer  and  produces  the  characteristic 
circular  fringe  pattern.  This,  however,  entails  a  loss  of  signal  power.  A 
plane  Fabry  Perot  basically  provides  a  transmission  that  depends  on  the  angle 
of  propagation  of  the  beam  with  respect  to  the  normal  to  the  reflectors.  It 
has  a  high  throughput  for  components  at  the  proper  frequency  and  angle.  The 
best  that  can  be  done  is  to  arrange  the  divergence  of  the  input  beam  so  that 
only  one  free  spectral  range  is  illuminated.  The  total  power  transmitted  in 
a  given  ring  is  then  given  by  the  total  power  incident  at  the  frequency  corre¬ 
sponding  to  the  ring  divided  by  the  finesse  of  the  interferometer. 

The  ease  with  which  a  high  resolving  power  car.  be  obtained  with  a  Fabry 
Perot  etalon  makes  it  attractive,  however,  and  it  will  be  used  for  a  demonstra¬ 
tion  of  the  frequency  domain  sampling  technique. 

The  most  convenient  form  of  spectrometer  for  use  in  the  system  ir  a 
grating.  Grating  can  have  a  high  throughput  and  produce  a  spatially  disfused 
out) ut  at  all  frequencies  simultaneously.  The  main  problem  is  that  of  ob¬ 
taining  sufficient  resolution.  Conventional  gratings  are  limited  to  a  resolu¬ 
tion  of  the  order  of  r00,000.  Echelle  gratings  are  capable  of  resolutions  in 
the  range  of  500,000  to  2,000,000.  For  any  type  of  multiple  beam  dispersive 
device,  the  resolution  will  be  determined  by  the  difierence  in  transit  time 
between  the  shortest  and  the  longest  path  in  the  device;  i,e., 
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for  a  grating  used  in  a  Littrow  arrangement  as  shown  in  Figure  It, 

2  W  sin  o 

At  - - 
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and 
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2  W 

-r-  sin  a 


Sche lie  gratings  are  typically  operated  with  a  diffraction  angle  a  ~  65°  so 
that  sin  m  0.9.  Thus  tc  achieve  a  resolution  of  106  at  a  wavelength  of 
1  micron,  a  grating  of  width 


2  sin  o 

is  required.  •  This  is  an  impractically  large  size  for  a  single  grating.  The 
required  resolution  may  be  obtained  by  multiple  diffractions.  Two- pass  and 
four- pass  instruments  have  been  fabricated  and  operated  successfully  For  a 
four-pass  instrument  the  grating  width  is  reduced  to  Ik  cm,  a  size  that  is 
commerically  obtainable  at  reasonable  cost. 

The  required  resolution  for  the  frequency  domain  sampling  system  is 
ob  ainable  easily  with  a  Fabry  Perot  filter  and  with  more  difficulty  with 
an  echelle  grating.  We  must  now  turn  to  the  problem  of  performing  the  required 
signal  processing  after  the  dispersion  is  obtained.  One  means  oi  performing 
s  would  be  to  implement  the  system  exactly  as  shown  in  Figure  7.  The  out¬ 
put  of  the  prating  could  be  imaged  upon  a  close  packed  array  of  detectors 
one  for  each  resoluable  frequency  clement.  The  output  of  each  detector  could 
ohen  be  modulated  in  an  electronic  mixer  with  the  appropriate  modulation  fre¬ 
quency  and  the  results  summed  electronically  to  yield  the  desired  output. 

There  is  no  fundamental  objection  to  this  approach  except  for  the  difficulty 
involved  in  maintaining  balance  among  all  the  separate  channels  and  the  large 
number  of  separate  electronic  channels  that  are  needed. 

Fortunately,  J he  entire  signal  processing  operation  after  dispersion  can 
be  conveniently  and  simply  accomplished  by  a  parallel  optical  processor.  One 
form  of  this  processor,  which  has  been  successfully  demonstrated,  is  tne  freouen-y 
multiplexed  Fabry  Perot  Interferometer  (Ref.  1).  This  device  simultaneously 
performs  the  detection,  modulation  and  signal  recombination  required  in  the 
present  case.  The  use  of  u  Fabry  Perot  interferometer  necessarily  entails  a 
loss  of  singal  as  discussed  above,  but  the  scheme,  may  be  extended  to  an  echelle 
grating  and  this  problem  can  be  eliminated. 
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The  form  of  the  Fabry  °erot  multiplexer  is  illustrated  schematically  in 
Figure  17a.  This  figure  illustrates  the  fringe  pattern  that  would  be  observed 
if  the  etalon  were  illuminated  by  a  line  spectrum  (in  this  case  only  four  lines 
are  shewn).  Figure  17b  shows  one  full  free  spectral  range.  To  use  this  output 
in  the  sampling  system  as  originally  described,  one  would  have  to  direct  each 
of  these  - ings  to  a  separate  detector  whose  output  is  subsequently  modulated. 
Schemes  to  do  this  involving  axicom  or  sets  of  annular  mirrors  have  been 
reported  in  Refs.  2,  3,  but  they  become  unweildy  for  a  large  number  of  channels. 
The  entire  problem  can  be  solved  by  modulating  the  rings  prior  to  detection. 

This  can  be  accomplished  by  the  so-called  multizone  disk  (Ref.  1).  The  disk 
is  shown  schematically  in  Figure  18.  It  consists  of  a  transparent  substrate 
that  is  divided  into  a  number  of  annular  zones.  The  mean  radius  of  each  zone 
is  centered  on  the  rings  that  are  observed  in  the  focal  plane  of  a  Fabry  Perot 
filter  that  is  illuminated  by  the  line  spectrum.  Each  zone  is  divided  into  an 
equal  number  of  transparent  and  opaque  area,  and  the  number  of  areas  increase 
as  the  radius  of  the  zone  increases .  Two  such  disks  are  placed  in  the  plane 
of  the  Fabry  Perot  fringes,  and  one  is  rotated  with  respect  to  the  other.  Each 
line  of  the  spectrum  is  thus  chopped  at  a  different  frequency  determined  by  the 
number  of  "teeth"  in  the  zone  and  the  rotational  rate.  The  entire  signal  that 
is  transmitted  through  the  device  is  focussed  onto  a  single  detector.  This  per¬ 
forms  an  operation  on  the  signal  that  is  very  nearly  the  one  required  to  recon¬ 
struct  the  signal.  Referring  to  Eq.  2.5,  ve  see  that  the  amplitude  of  each 
component  is  multiplied  by  a  function  sqr  (u^t  +  mCfet)  where  sqr  (u>t)  is  a 
square  wave  function  that  oscillates  between  0  and  1,  and  has  the  same  sign  and 
zero  crossings  as  cos  (cut).  The  derivation  leading  to  Eq.  2.9  follows  through 
exactly  as  before,  except  that  the  modulating  functions  are  square  waves  rather 
than  sin  waves.  If  we  arrange  the  disk  so  that  the  lowest  modulating  frequency 
is  greater  than  l/2  of  the  highest  modulating  frequency,  then  the  higher  harmonics 
of  the  square  waves  can  be  filtered  out  electrically  and  the  desired  output  will 
be  obtained. 

The  number  of  filter  elements  that  can  be  obtained  by  this  technique  is 
limited  only  by  the  finesse  of  the  interferometer.  The  maximum  modulation  fre¬ 
quency  is  determined  by  the  number  of  elements  that  can  be  placed  on  the  disk 
and  by  the  rotation  rate.  The  disk  can  be  made  by  drawing  it  on  a  large  scale 
with  a  computer  controlled  plotter  and  then  reducing  it  photographically.  It 
is  reasonable  to  assume  that  the  widths  of  the  teeth  in  the  zones  can  be  made 
as  small  as  20-30  microns.  For  an  inner  zone  of  1  cm  diameter,  this  would  give 
approximately  1000  elements.  If  the  disk  is  rotated  at  k)0  rps,  a  modulation 
frequency  of  400  KHz  can  be  obtained.  The  number  of  teeth  in  the  outermost 
zone  should  be  less  than  twice  the  number  in  the  inner  zone.  This  gives  a  total 
range  of  400  KHz.  If  there  are  lt-0  zones  on  the  disk,  the  frequencies  can  be 
spaced  by  10  KHz.  The  repetition  period  of  the  reconstructed  optical  signal 
coula  be  100  microseconds  in  this  case.  The  center  frequency  would  be  the  i.f. 
frequency  plus  or  minus  the  mean  modulation  frequency. 
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During  the  next  period  of  this  contract,  on  experimental  demonstration 
of  the  frequency  domain  sampling  using  the  multiplexed  Fabry  Perot  interfero¬ 
meter  -will  be  carried  out. 

A  preliminary  experiment  to  demonstrate  the  feasibility  of  individually 
heterodyning  the  line  components  cf  the  signal  and  the  reference  waver  has  been 
carried  out.  The  experimental  arrangement  was  similar  to  that  shown  in  Figure  5. 
A  scanning  Fabry  Perot  filter  was  placed  in  front  of  the  detector  so  that  the 
line3  of  the  spectrum  of  the  combined  signal  and  reference  could  be  examined 
individually.  The  detected  signal  at  the  i.f.  frequency  (40  MHz)  was  rectified 
and  displayed  on  an  oscilloscope.  T'  i  horizontal  axis  of  the  oscilloscope  was 
scanned  in  synchronism  with  the  spectrum  analyzer.  The  resulting  display  is 
shown  in  Figure  19a.  The  target  in  this  case  was  a  plane  mirror.  The  display 
shows  the  heterodyne  signals  from  each  of  the  lines  of  the  combined  signal  and 
reference  waveforms.  In  this  figure,  a  little  more  than  two  free  spectral  ranges 
are  shown.  The  oscillating  bandwidth  of  the  laser  was  somewhat  greater  than  the 
free  spectral  range  of  the  Fabry  Perot  (86  Hz)  so  that  there  is  some  overlapping 
of  orders.  When  an  additional  mirror  (~  65$  R)  was  placed  in  front  of  the  tar¬ 
get  mirror,  the  display  shown  in  Figure  19b  was  obtained.  The  two  mirrors  con¬ 
stitute  a  Fabry  Perot  filter.  When  operated  in  reflection,  the  filter  exhibits 
minima  at  frequencies  determined  by  Af  =  c/27  .  These  minima  may  be  seen  in 
the  figure.  The  location  of  the  minima  varied  with  the  spacing  of  the  mirrors 
as  would  be  expected  from  the  formula.  The  feasibility  of  heterodyning  the 
components  of  the ■ line  spectrum  individually  was  thus  confirmed. 

The  multiplexing  scheme  used  with  the  Fabry  Perot  can  be  applied  to  other 
types  of  dispersion  elements  such  as  an  echelle  grating.  A  possible  arrange¬ 
ment  is  shown  in  Figure  20.  The  output  of  the  grating  is  focussed  in  the  direc¬ 
tion  normal  to  the  dispersion  by  a  cylindrical  line.  This  line  is  then  scanned 
by  a  beam  deflector  across  a  mesh  whose  transparency  is  shown  schematically  in 
the  figure.  The  signal  transmitted  through  the  mask  is  then  focussed  into  a 
detector.  In  addition  to  eliminating  the  inefficiency  associated  with  the  Fabry 
Perot,  this  scheme  is  capable  of  modulating  at  much  higher  frequencies.  The 
modulation  frequency  is  limited  only  by  the  speed  of  the  beam  deflector  which 
may  be  mechanical,  electrooptical  or  acousto-optical. 
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